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Abstract. This paper is concerned with closed orbits of non-smooth 
vector fields on the plane. For a subclass of non-smooth vector fields we 
provide necessary and sufficient conditions for the existence of canard 
kind solutions. By means of a regularization we prove that the canard 
cycles are singular orbits of singular perturbation problems which are 
limit periodic sets of a sequence of limit cycles. Moreover, we generalize 
the Poincare Index for non-smooth vector fields. 



1. Introduction 

Piecewise-smooth systems are widespread within application areas such 
as engineering, economics, medicine, biology and ecology. The most com- 
mon piecewise-smooth systems involve either a discontinuity in the vector 
field, or in the orbit given by the integral solution x(t). In this paper we 
consider the former, that is, general systems where the vector field is in- 
dependently defined on either side of a smooth codimension one switching 
manifold. Three possible regions of the manifold are then apparent. At a 
crossing region the component of the vector field normal to the switching 
manifold has the same direction on both sides of the manifold (sometimes 
called sewing instead of crossing). At a stable sliding region both normal 
components of the vector field point toward the manifold. At an unsta- 
ble sliding region both normal components point away from the manifold. 
Piecewise-smooth systems with sliding are also known as Filippov systems. 
Clearly these three different scenarios lead to vastly different dynamics. An 
orbit that meets the switching manifold at a crossing region passes through 
it, but is non-differentiable at the crossing point. An orbit that impacts at 
a stable sliding region sticks becomes constrained (sticks) to the manifold. 
An orbit in an unstable sliding region slides along the switching manifold, 
but will depart it under any infinitesimal perturbation. Consequently, the 
only means by which a stable sliding orbit can escaping the switching man- 
ifold is tangentially, at the boundary of the sliding region. This leads to the 
observation that, under parameter variation, orbits in Filippov systems can 



1991 Mathematics Subject Classification. Primary 34C20, 34C26, 34D15, 34H05. 
Key words and phrases. Limit cycles, vector fields, singular perturbation, non-smooth 
vector fields, heteroclinic orbits, Poincare index, canard cycles. 



2 



C.A. BUZZI, T. DE CARVALHO, P.R. DA SILVA 



undergo a large variety of bifurcations, commonly called sliding bifurcations. 

In this paper we study piecewise-smooth system on open regions on the 
plane. Let IA C M. 2 be an open set and ECW given by £ = /~ 1 (0), where 
/ : IA — > K is a smooth function having £ 1 as a regular value (i.e. 
V/(p) 7^ 0, for any p G / -1 (0)). Clearly £ is the separating boundary of 
the regions £+ = {q G > 0} and £_ = {q € < 0}. We 

can assume that £ is represented, locally around a point q = (x,y), by the 
function f(x,y) = y. 

Designate by x r the space of C r vector fields on IA endowed with the 
C r -topology with r ) 1 or r = oo, large enough for our purposes. Call 
O r = £l r (U, f) the space of vector fields Xq :U\T, — > R 2 such that 



(1) X Q (x,y) 



Xt(x,y), for (x,y)€T, + , 
X 2 (x,y), for (x,2/)eE_, 



where Xj = (fi,gi) £ = 1,2. We write A"o = (Xi,X 2 ), which we will 
accept to be multivalued in the points of £. The trajectories of Xq are so- 
lutions of q = Xo(q), which has, in general, discontinuous righthand side. 
The basic results of differential equations, in this context, were stated by 
Filippov in [5]. Related theories can be found in [6j[9j 02] • I n this paper we 
consider finite discontinuities, i.e., the vector fields X\ and X 2 are defined in 
the set / _1 (0). Another kind of discontinuity of which the vector field tends 
to infinity when it approximates to the switching manifold can be obtained 
in the equations with impasse (more details in |10j). 

In what follows we will use the notation Xi.f(p) = (V/(p), Xi(p)). 
We distinguish the following regions on the discontinuity set E : 

(i) Si C E is the sewing region if (Xi.f)(X2-f) > on Si . 

(ii) £2 C E is the escaping region if (X\.f) > and {X 2 .f) < on £2. 

(iii) £3 C E is the sliding region if {X\.f) < and (X 2 .f) > on £3. 

Consider Xq € £l r . The sliding vector field associated to Xq is the vector 
field Xq tangent to £3 and defined at q € £3 by X^(q) = m — q with m 
being the point where the segment joining q + X\(q) and q + X 2 (q) is tan- 
gent to £3 (see Figure [T]). It is clear that if q € £3 then q € £2 for — Xq 
and then we can define the escaping vector field on £2 associated to Xq by 
Xq = — (— Xq) s . In what follows we use the notation Xq for both cases. 

Our main interest here is to study a special kind of typical minimal sets 
of non-smooth vector fields which will be called non-smooth " canard cycles 11 
(see Figure [2]) . A canard cycle is a graphic composed by pieces of orbit of 
X\, pieces of orbit of the sliding vector field Xq and/or pieces of orbit of 
X 2 . See Section 2 for a more precise definition. 
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FIGURE 1. Fillipov's convention. FIGURE 2. Canard cycle. 

An approximation of the non-smooth vector field Xq = (Xi, X2) by a 1- 
parameter family X t of smooth vector fields is called an e-regularization of 
Xq. We give the details about this process in section 4. A transition function 
is used to average X\ and X2 in order to get a family of smooth vector 
fields that approximates Xq. The main aim is to deduce certain dynamical 
properties of the non-smooth dynamical system from the regularized system. 
What is familiar may or may not be a matter of taste, at least it depends a 
lot on the dynamical properties of one's interest. The regularization process 
developed by Sotomayor and Teixeira produces a singular problem for which 
the discontinuous set is a center manifold. Via a blow up we establish a 
bridge between non-smooth systems and the geometric singular perturbation 
theory. 

Roughly speaking, the main results of this paper are the following: 

• In our first result (Theorem [1]) , for a subclass of non-smooth vector 
fields, we provide necessary and sufficient conditions for the existence 
of canard kind solutions. 

• In our second result (Theorem [2|), following the ideas exposed in 
[8], we prove that hyperbolic canard cycles are limit sets, according 
Hausdorf distance, of families of (smooth) hyperbolic limit cycles 
(this fact is not proved in [8]). The regularization process plus a 
blow up produce a singular perturbation problem V e . Our result 
implies that the canard cycle is the periodic limit set of closed or- 
bits of V e , with e — > 0. An open problem is to use the geometric 
singular perturbation theory proposed by Dumortier and Roussarie 
(center manifolds obtained via saturation by the flow plus blow up 
techniques, see [3] for details) to obtain the same result. 

• In our third result (Theorem [3]) we found an analogous for Poincare 
Index in the case of non-smooth vector fields. 

2. Preliminaries and statements of the main results 

Consider Xq £ Q r . We say that q € £ is a T,-regular point if 
(i) X 1 .f(q)X 2 .f(q) > or 
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(ii) X 1 .f{q)X 2 .f(q) < and A E (g) ^ (that is q £ £ 2 U ^3 and it is 
not a singular point of Xq). 

The points of £ which are not £-regular are called £- singular. We distin- 
guish two subsets in the set of S-singular points: £ c and £A We say that 
q £ £^ is a pseudo equilibrium of Xq if A E ((7) = and we say that (/ £ £ c 
is a ^-contact point if A^(q) 7^ and X\.f(q)X 2 .f(q) = (g is a contact 
point of Xq). 

A E-contact point q £ £ c is a T,-fold point of Xi if X\.f(q) = but 
Xff(q) 7^ 0. Moreover, g € £ is a visible (resp. invisible) T,-fold point of Xi 
if Xi./(g) = and Xf.f(q) > (resp. Xf.f(q) < 0). We say that q is a 
£-fold point of Xq if it is a £-fold point either of X\ or of X2. 

A pseudo equilibrium g £ £-^ is a T>-saddle provided one of the following 
condition is satisfied: (i) q £ £2 an d g is an attractor for A E or (ii) q £ S3 
and (; is a repeller for Aq\ A pseudo equilibrium q £ £ of Ao is a T,-repeller 
(resp. ^-attractor) provided q £ £2 (resp. g £ £3) and q is a repeller (resp. 
attractor) for Xq. A point q £ £ is a hyperbolic pseudo equilibrium of 
Ao if g is a hyperbolic equilibrium point of A E . 

Definition 1. Consider Xq £ Q r . 

(1) A curve T is a canard cycle ifT is closed and 

• r contains arcs of at least two of the vector fields Ai|s + , A2|e_ 
and A E or is composed by a single arc of X^; 

• the transition between arcs of X\ and arcs of X 2 happens in 
sewing points (and vice versa); 

• the transition between arcs of Ai ( or X 2 ) and arcs of Ag" hap- 
pens through T,-fold points or regular points in the escape or 
sliding arc, respecting the orientation. Moreover i/T ^ S then 
there exists at least one visible £-/oW point on each connected 
component of T n £ . 

(2) Let r be a canard cycle of Xq. We say that 

• r is a canard cycle of kind I if T meets £ just in sewing 
points; 

• T is a canard cycle of kind IlifT = £; 

• r is a canard cycle of kind III if T contains at least one 
visible T,-fold point of Xq. 

In Figures [7] an d appear canard cycles of kind I, II and III 
respectively. 

(3) Let r be a canard cycle. We say that T is hyperbolic if 

• r is of kind I and rj'(p) ^ 1 where r\ is the first return map 
defined on a segment T with p £ T fh 7; 

• r is of kind II; 

• r is of kind III and or T D £ C £x U £2 or T n £ C £^ U £3. 
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FIGURE 3. Canard 
cycle of kind I. 



FIGURE 4. Canard 
cycle of kind II. 



In p] is proved that the e-regularization of non-smooth vector fields Xq 
with hyperbolic canard cycles has hyperbolic limit cycles. 

Definition 2. Let {A B)x t be an arc of Xi joining the visible E -fold point A 
to the point B = Xi iti E. We say that {AJ^)x i has focal kind if there is not 
Yi-fold points between A and B (see Figure^) and we say that (A~i)xi 

has 

graphic kind if it has only one E-/oW point between A and B (see Figure 
W, * = 1,2. 





A B 

FIGURE 5. Focal kind arc. 



A B 

FIGURE 6. Graphic kind arc. 



We remark that for a hyperbolic canard cycle we have that each connected 
component of V H E has only one E-fold point (See [8] for more details). 
By using the previous notation, our results are: 

Theorem 1. Let Xq = (Xi,X 2 ) € £l r be a non-smooth vector field with Xq 
presenting only one E -fold point A which is visible. Denote 71 the arc of Xi 
(i = 1 or i = 2) which passes through A and call B the transversal contact 
point 0/71 with E. Then Xq has a canard cycle T if and only if the following 
conditions are satisfied: (i) the component 71 of T which passes through A is 
a focal kind arc; (ii) X1/.X2/ < in (A, B] and (Hi) {X±,X2} is a linearly 
independent set in [A, B] . Moreover, T is of kind III. 

Theorem 2. Let Tq be a hyperbolic canard cycle of Xq. Then for any 
e > the regularized vector field X e , has a hyperbolic limit cycle T e such 
that T e — >• Tq when e — > 0. 

We remark that the Hausdorff distance between compact sets of M? is: 
D(K X ,K 2 )= max {d(z 1 ,K 2 ),d(z 2 ,K 1 )}. 
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Theorem 3. Let Tq be a hyperbolic canard cycle of the non-smooth vector 
field Xq. If {p\, . . ., pk} is the set of fixed or pseudo equilibrium points (all 
hyperbolic ) of Xq inside Tq then the index of Tq with respect to Xq is the 
sum of the index of pi, for i = 1, k. Moreover, this sum is equal to one. 

In section 5 we will define index of non-smooth vector fields. 

The paper is organized as follows. In Sections 3, 4 and 5 we prove The- 
orems [H [2] and El respectively. In section 6 we apply Theorem [T] to study 
a class of non-smooth vector fields Xq € f2 r with just one focal kind arc 
and its bifurcation and we use the singular perturbation theory to study 
hyperbolic canard cycles. 



3. Proof of the Theorem CD 
In this section we prove the first result of the paper. 

Proof. First we prove that (i),(ii) and (iii) imply the existence of the canard 
cycle. Since X\f .X<if < in (A, B] the piece of £ between A and B is part 
of a escaping region or a sliding region. Moreover since {X%, X2} is a linearly 
independent set in [A, B] the system does not have pseudo equilibrium points 
in [A, B\. Without lost of generality, [A, B] is part of the sliding region like 
in Figure [3 The curve T = ^/i L) [B , A] is a hyperbolic canard cycle of kind 
III. We remark that this canard cycle takes place in just one side of £. 

Now we prove that (i),(ii) and (iii) are necessaries conditions for the exis- 
tence of this particular kind of canard cycle. Since T is a hyperbolic canard 
cycle of kind III with just one E-fold point, T takes place in just one side 
of S. In fact, if it does not occur, then T returns to E at least twice and 
so there exists at least a second S-fold point. Without lost of generality we 
suppose that T is on the side corresponding to X\. We denote by 71 the part 
of the cycle T which is a trajectory of X\. Thus we have that 71 is a focal 
kind arc because if it is a graphic kind arc then there is another S-fold point 
on (^4, B) (see Figure [7]). Since V has no one arc of X2, the point B belongs 
to an escaping region or a sliding region and so X\f{B).X2f{B) < 0. Let 
us assume that B E S3. Since 71 meets £ in the point B, the flow slides via 
X^ 1 until the point A because there are not another S-fold point between A 
and B; therefore X\f.X2f < in (A, B\. Moreover, the linear independence 
of {Xi,^} on [A, B] follows from the non-existence of pseudo equilibrium 
points on [A, B\. 

Now, we will define an auxiliar function which will be useful in the sequel. 

Take (^4, B) C S2 U S3 contained in the escaping or in the sliding region. 
In (A, B) consider the point C = (61,62), the vectors X\(C) = {Di^D^) 
and X2{C) = (E±,E2) (as illustrated in Figured]). The straight segment 
passing through C + X\(C) and C + -X^C) meets S in a point p{C). We 
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71 



x 2 

FIGURE 7. Canard cycle with just one E-fold point, 
define the C r -application 

p: (A, B) — > E 
z i — > p(z). 

We can choose local coordinates such that E is the x-axis; so C = (C\,0) 
and p{C) elx {0}. The direction function on E is defined by 

H : (A, B) — ► R 

z i — > p{z) — z. 

Xi sjC + X 2 {C) 

s — i ^—Mm, 

A C\/ 5 

X 2 C + Xi(C) 

FIGURE 8. Direction function. 

We have that is a C r -application and 

• if H(C) < then the orientation of Xq in a small neighborhood of 
C is from B to ^4; 

• if H{C) = then C € S^; 

• if H(C) > then the orientation of Xq in a small neighborhood of 
C is from A to -B. 

Simple calculations shows that p(d) = {Dl+Cl){ ^ ) Z[^ ) {El+Cl) ■ 

Assuming all the hypothesis of Theorem Q] we have the following corollary. 

Corollary 1. The non- smooth vector field Xq has a canard cycle V if and 
only if the direction function H : [A, B] — > K. is a well defined function and 
it has no zeros. Moreover, T is of kind III. 

4. Proof of Theorem [2] 

First of all we present the concept of e-regularization of non-smooth vector 
fields. It was introduced by Sotomayor and Teixeira in [9]. The regulariza- 
tion gives the mathematical tool to study the stability of these systems, 
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according to the program introduced by Peixoto. The method consists in 
the analysis of the regularized vector field which is a smooth approximation 
of the non-smooth vector field. Using this process we get a f-parameter 
family of vector fields X e € x r (-^>^ 2 ) such that for each eo > fixed we 
have 

(i) X eo is equal to X\ in all points of S + whose distance to £ is bigger 



(ii) X eo is equal to X2 in all points of S_ whose distance to £ is bigger 



Definition 3. A C°° function ip : R — > R is a transition function ifip(x) = 
— 1 for x ^ —1, ip(x) = 1 for x 1 and f'(x) > if x 6 (— 1, 1). The e- 
regularization of Xq = (X\,X2) is the 1-parameter family X e G x r given 
by 



with ip e {x) = ip(x/e), for e > 0. 

In order to prove Theorem [2] we need to construct a special neighborhood 
of arbitrary diameter for hyperbolic canard cycles. 

Construction of a neighborhood of diameter fi around a hyperbolic 
canard cycle. Here we describe a method to construct a tubular neighbor- 
hood of diameter [i around a hyperbolic canard cycle. This presentation is 
done for canard cycles of kind III, but the ideas can also be extended for 
kinds I or II. We will be particularly interested in two of them: the ones 
that take place on just one side of £ and with just one visible S-fold point 
and the ones that take place on the two sides of £ with two visible S-fold 
points (one for X\ and another one for X2). 

Case 1- One T,-fold point. Denote by T the hyperbolic canard cycle of 
kind III with just one S-fold point and with orientation showed in Figured] 
(the reverse orientation is treated in a similar way). Consider the strip of 
diameter \x around £ = {y = 0}. Let p\ and q\ be points in {y = fi} n T. 
Take an arc 72 of the vector field X\ passing to the point P2 € {y = /i} in 
such a way that P2 stays on the left of p\ and such that 72 returns to the line 
y = n in a point q2 which is in a neighborhood of q\. Take this trajectory 
satisfying d(T, 72) < ^ (this is possible by the continuity of X\). Let r\ be 
the point where the arc of X\ through by p\ first meets the straight line 
y = [i for negative time. Analogously take an arc 71 of the field X\ passing 
by the point T2 in such a way that T2 stays on the left of r\ and such that 
71 has second return to y = fi in a point q$ which is in a neighborhood of 
q\. Take this trajectory satisfying <i(r,7i) < ^. On S, on the left of the 
S-fold point A, the flow of X\ is oriented to up, so it is possible to construct 
a transversal section 02 joining P2 to the straight line y = in such a way 
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that the same trajectories of X\ cross transversally o~i and the segment 
pypi- Take ti the point where 01 meets the straight line y = satisfying 
d(T, a 2) < /U, as before. Moreover, on E, on the right of the point B, the flow 
of Xi is oriented for down, so it is possible to construct a transversal section 
62 joining 52 to the straight line y = in such a way that the trajectories 
of X\ that cross transversally 62 do not cross the segment q± q2- Let S2 be 
the point where O2 meets the straight line y = (here we also need to take 
care for d(r,# 2 ) < Since [A, B] is a sliding region, the flow of X 2 is 
transversal to E. In the straight line y = —/j,, consider the points U2 and 
i>2 (v,2 is on the left of ^2) satisfying that the trajectories of X2 crossing the 
transversal sections A2 = ^2 ^2 and 62 = S2 v 2 meet transversally E at the 
segment £2 ^2 (again, we need to take care for d(T, A2) < p and d(T, 62) < 




FIGURE 9. Tubular neighborhood of a canard cycle with one S-fold. 

In this way, the strip defined by the closed curve 71 U 93^2 and by the 
closed curve 72 U 82 U 62 U v 2 ^2 U A2 U 02 is a tubular neighborhood of T of 
diameter [i. Note that the flow of Xq = (Xi,X2) is arriving in this neigh- 
borhood and never it departs from it. 

Case 2- Two Yj-fold points. Now we study the hyperbolic canard cycles 
of kind III with two visible E-fold points, being one for X\ and the other 
one for X2 , like showed in Figure [TTJJ We work with canard cycles T that 
have only escaping regions on E = {y = 0} (the case with sliding regions 
is treated similarly). Consider the strip of diameter fi around E. Let p± 
and q\ be points in {y = /j,} n T. Take an arc 71 of the vector field X\ 
through t\ € {y = /j,} satisfying that t\ stays on the left of pi and such that 
71 returns to the line y = fi in a point u\ which is in a neighborhood of 
q\. Take this trajectory satisfying that d(r,7i) < ^. Take an arc a\ of the 
vector field X\ through v\ 6 {y = /i} satisfying that v\ stays on the right of 
pi and such that o~\ has second return on the straight line y = p, in a point 
x±, even take this trajectory with the particularity that d(T,ai) < ^. We 
repeat the same argument for the vector field X2 and we found the points 
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P2, Q2, t2, U2, V2 and X2 respectively, and the curves 72 and 02. Let c be the 
point on En 02, d be the point on S D 71, e be the point on S D 72 and / 
be the point on £ n o\ as indicated in Figure [TOl On E, take the points g 
on the left of c, h between A and d, i between e and A' and j on the right 
of /; satisfying that the arcs 6\ (joining g to 27), p\ (joining u\ to h), rj\ 
(joining i to t{) and 7Ti (joining v\ to j) are transversal sections for X\ and 
the arcs 62 (joining g to V2), P2 (joining h to t^), 7/2 (joining U2 to i) and 7r2 
(joining j and X2) are transversal sections for Xq with the distance from V 
to any one of this arcs less than p. 



FIGURE 10. Tubular neighborhood for canard cycles of kind III 
with two S-fold points. 

In this way, the strip defined by the closed curve 71 U p\ U P2 U 72 U 772 U rj\ 
and by the closed curve a\ U 6\ U 62 U 02 U 7r2 U 7Ti is a tubular neighborhood 
for r of diameter p. Note that the flow of Xq = (X\,X2) is departing from 
the tubular neighborhood and never it arrives in it. 

• Since this neighborhood bounds a region where the non-smooth vec- 
tor field Xq = (Xi,X2) is arriving in or it is departing from them, 
it makes sense to say attractor canard cycle or repeller canard cycle. 

• In the neighborhoods constructed before we allow that trajectories 
can make part of them, however it is possible to do it with the flow of 
Xq being transversal to the boundaries of the tubular neighborhoods. 
In fact, it is enough to replace the trajectories by transversal curves. 
It is important for the construction of the tubular neighborhood of 
the canard cycles of kind I. Thus we make a construction like we 
made before but now we can use for this, the first return application 
n and thus if 7/ < 1 we have an attractor canard cycle and if 7/ > 1 
we have a repeller canard cycle. 

• For canard cycles of kind II is enough to take the strip of diameter 
p in the beginning of the construction as the tubular neighborhood. 
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• Any other hyperbolic canard cycle is an arrangement of pieces of 
the canard cycles described above and so we can construct a tubular 
neighborhood for it arranging the previous tubular neighborhoods. 

Proof of Theorem [2j Let Fq be a canard cycle of Xq and let V e be a 
tubular neighborhood of diameter e around Tq. Since Xq is transversal to 
the boundary of V e , by continuity, the regularized vector field X e also is 
transversal to the boundary of V t . Assume that Tq is an attractor canard 
cycle, so the flow of Xq is arriving in the neighborhood V e and consequently 
the flow of X e also is arriving in the neighborhood V t . As there are not 
fixed points in V e , applying the Poincare-Bendixson Theorem we conclude 
that there exists an attractor limit cycle T e inside V t . Moreover with a more 
detailed analysis we can prove that it is hyperbolic (see [8] for instance). 
Since every paths that compose V e depends continuously of e we have that 
the diameter of the tubular neighborhood is a continuous function of the 
variable e. Therefore making e — > we conclude that T e — >• Tq (see Figure 

HB. | 




FIGURE 11. Cycles convergence. 

We remark that if Tq is an attractor (resp. repeller) hyperbolic canard 
cycle of Xq, then the same occurs for T e and X e . 

5. Proof of Theorem [3] 

Now we start our discussion about the third result. Let / = [0, 1] be an 
interval and a : I U be an oriented closed continuous path. Suppose 
that there are no critical points of Xq on a. Let us move a point P along 
the curve in the counterclockwise direction. The vector Xq(P) will rotate 
during the motion. When P returns to its starting place after one revolution 
along the curve a, Xq(P) also returns to its original position. During the 
journey Xq(P) will make some whole number of revolutions. Counting these 
revolutions positively if they are counterclockwise, negatively if they are 
clockwise, the resulting algebraic sum of the number of revolutions is called 
the index of a with respect to Xq, and is denoted by I{Xq,o). 
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To calculate I(Xo,o~) is convenient normalize X(a(t)) as an unit vector 
at the origin. In this way, we can define a function 9 : I — )■ R such that 

- lim iiyi^SIn =- lim (cos0(t),sin0(t)) 
t-^-t- \\X(a(tj)\\ t^rt- 

for every t G I. The function 9 is called angle function. 

We observe that in the case of smooth vector fields, the angle function 

is always continuous, but in the case of non-smooth vector field it admits 

a "jump" when the path pass to a point Sj G X, i G N. Therefore, we 

establish a rule for this jump; at Si = o~(ti) the angle function oscillates 

from lim 9{t) to lim 9{t). If lim 9(t) G I k -i = (2(k - l)7r,2fevr) and 
t-+tr t->tf t^t~ 

lim 9(t) G Ik = (2kir, 2(k + l)vr), where k € Z, we add 1 to the the number 
I(Xo,<t); if lim G and lim G Ik-i we add —1 to the number 
I(Xo,o~). We always consider that the jump of the vector lim X(a(t)) to 
the vector lim X(a(t)) occurs by the smallest angle between this vectors. 




FIGURE 12. Angle function: at si it has a jump of size ct\ and 
at S2 it has a jump of size ot2 ■ 

The difference 9(1) — 9(0) is a multiple of 27r, and 

T( X „1 g(l)-g(0) 
/(X ,a) = 

is an integer independent of the chosen cr-parametrization. This number 
also is called Poincare Index of the curve a with relation to the non-smooth 
vector field Xq. 

Our interest here is to calculate the index of canard cycles surrounding 
fixed or pseudo equilibrium points that are the critical points of Xq. We will 
see that, different from the smooth case, given two canard cycles T\ and T2 
surrounding the same critical points of Xq we have I(Xq, T±) ^ I(Xq,T2) in 
general. 



Example. Consider the configuration described in Figure [T3j Let 71 be an 
arc of X\ joining the S-fold points a and b; e the S-fold point of X2; c and 
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d points in the escaping region eb; g, h and i points in the sliding region 
fa where / is a invisible X-fold point of X\\ 62 the arc of X2 joining c and 
h, Si the arc of X\ joining d and g; and 72 the arc of X2 joining e and i. 
Consider the fixed point of X2 named by P2 in Figure [131 We can choose 
closed curves such that the Poincare index to non-smooth vector fields is any 
natural number. For example, take the path /iaU7iU?>cU<52 and the index 
is 1; in an analogous way if we take the path iaU7iU6cU<52U/iaU7iU6eU72 
the index is 2. Repeating this argument we can found closed curves such 
that the Poincare index is any integer number that we wish. 




FIGURE 13. Non- hyperbolic canard cycle. 



The previous situation just occurs because the canard cycle given, for 
example, by 71 U be U 82 U ha, is non-hyperbolic because in its composition 
we can found pieces of escaping region and pieces of sliding region. If we 
eliminate this possibility we have the next theorem, which has an analogous 
in the case of smooth vector fields. 

Remark 1. We recall that if Xq has a hyperbolic S-saddle (or a hyperbolic 
S-focus) so then the regularized vector field X t has a hyperbolic saddle (hy- 
perbolic focus) s e where s e — > sq when e — > (for details see [8]). Moreover, 
we can verify that if sq is a saddle or a S-saddle then we can take a suffi- 
ciently small closed path a around sq and prove that I{Xq, a) = — 1 (and if 
so is a focus or a S-attractor or a S-repeller then we can take a sufficiently 
small closed path a around so and prove that I(Xq, a) = 1). When the path 
a is sufficiently small to have just one critical point of Xq, named so, in its 
interior we use the notation I So (Xo,a) to denote its Poincare Index. 

5.1. Proof of Theorem [3l We want to prove that if To is a hyperbolic 
canard cycle of the non-smooth vector field Xq and if p\ , ■ ■ ■, Pk are the only 
ones critical points (all hyperbolic) of Xq inside Tq then I(Xq,Tq) is well 
defined and 

k 

I(X ,T ) = ^/ Pi (x ,r ) = 1. 

i=l 
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First of all we assume that the index is well defined. Let Xq be a non-smooth 
vector field with a hyperbolic canard cycle Tq and k hyperbolic critical points 
of Xq inside Tq. Thus, the regularized vector field X e has a hyperbolic limit 
cycle T e and k hyperbolic fixed points inside T t . So, by the Poincare Index 
Theorem for smooth vector fields the index calculated in T e in relation to 
X e is the sum of the index of the fixed points of X e inside r e and this sum 
is equal to 1. Since T e — > Tq we conclude that the index calculated in Tq in 
relation to Xq is the sum of the index of the critical points of Xq inside Tq 
and this sum is equal to 1 (see remark [IJ) . In order to finish the proof we 
verify that the index is well defined in the case that the closed curve Tq is a 
hyperbolic canard cycle. In fact, let Tq be a hyperbolic canard cycle of Xq. 
Let us assume that Tq is of one kind described in Figures [HI or [15] below. 




FIGURE 14. p x is repeller. FIGURE 15. pi is attractor. 

We will prove that there is not danger of ambiguity in the choose of the 
closed paths, differently that what happens in the previous example. In 
Figure [T4l we consider the hyperbolic canard cycle given by Tq = 71 U B A. 
Any canard cycle of Xq having pieces of sliding region must to pass by the 
E-fold point A, now walk in 71 (which is the only one possibility that we 
have!) and meet the point B on E. The unique choice we have is return to 
A, closing the path, without ambiguity. For Figure [151 the analysis is more 
interesting. Obviously we can use the trick of take the vector field — Xq 
and obtain an analogous result that the previous, however we prefer give 
here a complete idea to the case described in Figure [15j Since in the semi 
straight line r = A% we have an escaping region, the canard cycle must have 
only escaping region in its composition. Note that, if we choose to depart 
from r by a point in the segment (A, B) then this path does not return to 
(A,B) C E (it will move spirally around the focus pi), if we choose going 
out from r by a point after B in r then this path also does not return to 
(A, B) C E (to this path return to (A, B) it must return in a sliding region, 
what is not allowed because only escaping regions compose this hyperbolic 
canard cycle). So we must leave r by the point B and to close the curve. 
Therefore in any case there is not danger of ambiguity in the choose of the 
closed curve and so the Poincare index for non-smooth vector fields is well 
defined. To hyperbolic canard cycles of kind III with another particularities 
is enough repeat the ideas exposed here. To hyperbolic canard cycles of 
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kinds I and II clearly there is not danger of ambiguity in the choose of the 
closed curves once there are not escaping or sliding region in its composition. 

I 

Corollary 2. Under the hypothesis of the previous theorem and assuming 
that all canard cycles of Xq are hyperbolic, we have that: 

(1) If Tq is a canard cycle then inside Tq there exist (2n + 1) critical 
points of Xq, being n saddles or T,-saddles and (n + 1) focus, S- 
repeller or T,-attractor. 

(2) If all critical points of Xq are saddle or Y,-saddle then Xq does not 
have canard cycles. 

Proof. Since the index of each saddle and each E-saddle point is equal to 
— 1 and the index of each other critical point of Xq is equal to 1 the result 
is an immediate consequence of Theorem [3J 

6. Applications and Examples 

6.1. Heteroclinic Orbits. Consider the notation of the Theorem [TJ We 
give now an example of a curve that satisfies all the hypothesis in this theo- 
rem except that there exists a point C € (^4, B) such that the vectors X\(C) 
and X2(C) are not linearly independent; instead of V obtained in the theo- 
rem we have here a "T,-loop", that is a S- saddle- attractor with connection 
between E-separatrices. 

Example. Consider the non-smooth vector field Xq = (Xi, X2) 
with Xi(x, z) = (x + z — 1, —x + 2 — X2(x, z) = (—x 2 + |x — |, 1) and 
discontinuity set given by the x-axis, i.e., f(x,z) = z. On z = 0, we have 
Xi(x, 0) = (x — 1, —x — 1) and X 2 (x, 0) = {—x 2 + |x — |, 1) and so, 

(3) (X 1 .f)(x,0) = -x-l, (X 2 .f)(x,0) = l. 

In this way, we can conclude that x = — 1 is a S-fold point of X\ which 
determines a focal kind arc. For x > — 1 we have that S is a sliding region 
and for x < —1 it is a sewing region (see Figure ITH]) . We show now that there 
exists a point C in the semi straight line x > — 1 for which X\(C) = \.X2(C). 
In fact, ifXi(x,0) = X.X 2 (x,0) then h(x) = (x-l) 2 (x+|) = 0. The graphic 
of h{x) is given in Figure [171 We observe that h is equal to — H, where H is 
the direction function defined previously. So, we have the situation described 
in the Theorem [H except that X\ and X2 are not linearly independent in 
x = 1 where X^ 1 has an equilibrium point. The orientation of X^ is in 
direction to the E-fold point because for x = | we have X\(^, 0) = (^-, ^), 
^2(5,0) = (0, 1) and so, the direction function H is negative (H(^) = — |), 
analogously for x = | we have Xi(|,0) = (^,^), ^2(^,0) = (^-,1) and 
so, the direction function H also is negative (-f^(f) = — j^)- The pseudo 
equilibrium p = (1,0) is a S-saddle-attractor where the S-separatrices are 
connected. 
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FIGURE 16. £-saddlc-attractor 



Figure 

17. Graphic of 
h. 



with E-separatrices connection. h. 

Example (Bifurcation of the previous example) In the previous example, 
the pseudo equilibrium with a "loop" was found because the function h has 
a double-zero at x = 1. So, we can conclude that putting small variations on 
the fields X\ and X2 a new function h appears, with two simple real zeros 
(see Figure I2T1) or without real zeros (see Figure fT9j) in a neighborhood of 
x = 1. The phase portrait of Xo for small variations on X\ and X2 are 
showed in Figures [18] and [20j Note that for the case showed in Figure [T8| 
the direction function H is always negative in a neighborhood of x = 1 
and we can apply the corollary Q] to conclude that there exists a hyperbolic 
canard cycle of kind III and, for the case showed in Figure [20j the direction 
function H has two simple real zeros in a neighborhood of x = 1 and assumes 
positive values between this two points and negative values in the rest of 
the sliding region. In this way we have a bifurcation model where imposing 
small variations we can have a hyperbolic canard cycle or we can have two 
pseudo equilibrium points with a stable connection between its separatrices. 
We call this a S-Loop Bifurcation. 



Canard Cycle of kind III. 

Following the notation of Theorem Q] and the ideas just exposed we can 
state the next Proposition: 

Proposition 1. Let Xq a non-smooth vector field. 




FIGURE 18. Hyperbolic 



FIGURE 19. Graphic 
ofh. 
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Figure 20. E-saddlc Figure 21. Graphic 

and E-attractor with 
E-separatrices connection. 

(1) Xq has an unstable configuration topologically equivalent to that one 
in Figure [Td\ if and only if (i) the direction function H is well defined 
in [A,B], (ii) H has a single zero in (A,B) and (Hi) H(B) < 0. 

(2) Xq has a stable configuration topologically equivalent to that one in 
Figure [2B if and only if (i) the direction function H is well defined 
in [A, B], (ii) H has a single zero in (A,B) and (Hi) H(B) > 0. 

(3) Xq has an unstable configuration topologically equivalent to that one 
in Figure {23\ if and only if (i) the direction function H is well defined 
in [A,B\, (ii) H do not have zeros in (A,B) and (Hi) H(B) = 0. 

Proof. It is straightforward following what is done in the previous example. 




FIGURE 22. Stable configuration. FIGURE 23. Unstable configuration. 

Moreover, concerning with the item (1) of the Proposition 1 small per- 
turbations in Xq produces the effects showed in the previous example and 
concerning with the item (3) of the Proposition 1 small perturbations in Xq 
produces effects described in Theorem [H in the previous example and in the 
item (2) of the Proposition 1. 

6.2. Canard Cycles and Singular Perturbations Problems. In this 
section we show how the regularization process gives a singular perturba- 
tion problem. In this context, the canard cycles defined in this paper can be 
considered as limit periodic sets of singular problems. First of all we present 
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some basic definitions. 

Definition 4. Let U C R 2 6e an open subset and take e ^ 0. A singular 
perturbation problem in U (SP-Problem) is a differential system which can 
be written like 

(4) x' = dx/dr = l(x,y, e), y' = dy /dr = em(x,y, e) 
or equivalently, after the time re-scaling t = er 

(5) ex = edx/dt = l(x,y,e), y = dy/dt = m(x,y,e), 
with (x,y) € U and l,m smooth in all variables. 

The understanding of the phase portrait of the vector field associated 
to a SP-problem is the main goal of the geometric singular perturbation 
theory (GSP-theory). The techniques of GSP-theory can be used to obtain 
information on the dynamics of @) for small values of e > 0, mainly in 
searching limit cycles. System is called the fast system, and © the slow 
system of SP-problem. Observe that for e > the phase portraits of the 
fast and the slow systems coincide. For e = 0, let <S be the set 

(6) S = {(x,y):f(x,y,0) = 0} 

of all singular points of ([!]). We call S the slow manifold of the singular 
perturbation problem and it is important to notice that equation © defines 
a dynamical system, on 5, called the reduced problem: 

(7) f(x,y,0) = 0, y = g(x,y,0). 

Combining results on the dynamics of these two limiting problems, with 
e = 0, one obtains information on the dynamics of X e for small values of e. 
We refer to [1] for an introduction to the general theory of singular pertur- 
bations. Related problems can be seen in [1], [3] and [11]. Let us apply the 
techniques of GSP-Theory to study hyperbolic canard cycles. 

Example. Consider the non-smooth vector field Xq = (Xi,X2) with 
Xi(x,y) = (x + y - 1,-x + y + 1), X 2 {x,y) = (1,2) and f{x,y) = x. 
The regularized vector field becomes 

x = ^ + ^(f)^, 

y = ^±3 + ^(f)^±ti. 

where ^(f ) is the transition function. Making the change of variables x = 
r. cos 9 and e = r. sin 9 we obtain 



sing ( j- cos 2 e +y +ip (cot flp 



cos 9+y—2 
2 



(8) 

y = ~ r - cos °+y+ 3 + y (cot fl) ~ r - cob g+y-i . 

In the blowing up locus r = the fast dynamics is determined by the 
system 



sin^l + ^cotfl)^-^) , y' = 0; 
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and the slow dynamics on the slow manifold is determined by the reduced 
system 

! + ^(cot 0)^ = 0, y = ^ + ^(cot0)^. 

We remark that the slow manifold is implicitly defined by |+(^(cot 6)^-^- = 

and y(6) defined in this way is such that lim y(9) = 1, lim y(9) = — oo. 

9 — >\ e — 

In the phase portrait on the blowing up locus double arrow over one the 

trajectory means that the trajectory is of the fast dynamical system, and 

simple arrow means that the trajectory is of the slow dynamical system. So, 

we can draw the slow variety and its orientation and give the orientation of 

the fast flow (see Figure |24"|) • 




FIGURE 24. Singular FIGURE 25. Singular 
perturbation of a Ca- perturbation of a Ca- 
nard with one S-fold nard with tw0 s - fold 
point. points. 



In the final example of [2] the authors apply the GSP-Theory to the non- 
smooth vector field X (x, y) = (X 1 (x, y), X 2 (x, y)) = ((3y 2 -y-2, 1), (-3y 2 - 
y + 2, —1)) and obtain the SP-problem which behavior is described in Figure 
[25J 

In (7j the authors prove that S 2 n S 3 is homeomorphic to the slow variety 
and that the sliding vector field X^ is topologically equivalent to the reduced 
problem. So, we can apply step-by-step the method described in section 4 of 
this paper and found tubular neighborhoods for the canard cycles in Figures 
l24l and l25l Moreover, we can apply the Theorem [2] and conclude that it is a 
limit set (making r —¥ 0) of hyperbolic limit cycles. This also is true to any 
one hyperbolic canard cycle. 

Acknowledgments. The first and the third authors are partially sup- 
ported by CAPES and CNPq, the second author is partially supported by 
a FAPESP-BRAZIL grant 2007/08707-5. 



20 



C.A. BUZZI, T. DE CARVALHO, P.R. DA SILVA 



References 

[1] C.A. Buzzi, P.R. DA Silva AND M.A. Teixeira, Singular Perturbation Problems 
For Time Reversible Systems, Proc. Amer. Math. Soc, 133 (2005), 3323-3331. 

[2] C.A. Buzzi, P.R. da Silva and M. A. Teixeira, A Singular Approach To Discon- 
tinuous Vector Fields On The Plane, Journal of Differential Equations, 231 (2006), 
633-655. 

[3] F. Dumortier AND R. Roussarie, Canard Cycles And Center Manifolds, Memoirs 
Amer. Mat. Soc. 121, 1996. 

[4] N. Fenichel, Geometric Singular Perturbation Theory For Ordinary Differential 
Equations, Journal of Differential Equations 31 (1979), 53-98. 

[5] A.F. Filippov, Differential Equations With Discontinuous Righthand Sides, Mathe- 
matics and its Applications (Soviet Series), Kluwer Academic Publishers-Dordrecht, 
1988. 

[6] V. S. Kozlova, Roughness Of A Discontinuous System, Vestinik Moskovskogo Uni- 

versiteta, Matematika 5 (1984), 16-20. 
[7] J. Llibre, P.R. Silva and M.A. Teixeira, Sliding Vector Fields Via Slow-Fast 

Systems, Bulletin of the Belgian Mathematical Society Simon Stevin 15-5 (2008), 

851-869. 

[8] J. Sotomayor and A.L. Machado, Structurally Stable Discontinuous Vector Fields 

On The Plane, Qual. Theory of Dynamical Systems, 3 (2002), 227-250. 
[9] J. Sotomayor and M.A. Teixeira, Regularization Of Discontinuous Vector Fields, 
International Conference on Differential Equations, Lisboa (1996), 207-223. 

[10] J. Sotomayor and M. Zhitomirskii, Impasse Singularities of Differential Systems 
of the Form A(x)x' = F{x), J. Diff. Equations 169, n°2, (2001), 567-587. 

[11] P. Szmolyan, Transversal Heteroclmic And Homoclmic Orbits In Singular Pertur- 
bation Problems, Journal of Differential Equations 92 (1991), 252-281. 

[12] M.A. Teixeira, Generic Singularities Of Discontinuous Vector Fields, An. Ac. Bras. 
Cienc. 53, n°2, (1991), 257-260. 

1 2 3 IBILCE-UNESP, CEP 15054-000 S. J. Rio Preto, Sao Paulo, Brazil 
E-mail address: 1 buzzi@ibilce .unesp .br 2 ti-car@hotmail . com 
E-mail address: 3 prs@ibilce .unesp .br 



